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ABSTRACT 



The asymptotic giant branch (AGB) stars X Her and TX Psc have been imaged at 70 and 160 /im with the PACS instrument onboard 
the Herschel satellite, as part of the large MESS (Mass loss of Evolved StarS) Guaranteed Time Key Program. The images reveal an 
axisymmetric extended structure with its axis oriented along the space motion of the stars. This extended structure is very likely to 
be shaped by the interaction of the wind ejected by the AGB star with the surrounding interstellar medium (ISM). As predicted by 
numerical simulations, the detailed structure of the wind-ISM interface depends upon the relative velocity between star+wind and the 
ISM, which is large for these two stars (108 and 55 km s"' for X Her and TX Psc, respectively). In both cases, there is a compact blob 
upstream whose origin is not fully elucidated, but that could be the signature of some instability in the wind-ISM shock. Deconvolved 
images of X Her and TX Psc reveal several discrete structures along the outermost filaments, which could be Kelvin-Helmholtz 
vortices. Finally, TX Psc is surrounded by an almost circular ring (the signature of the termination shock?) that contrasts with the 
outer, more structured filaments. A similar inner circular structure seems to be present in X Her as well, albeit less clearly. 
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1. Introduction 

The MES S (Mass loss of Evolved StarS) Guaranteed Time Key 
Program dGroenewe^en et al.l201 1) observed with the PACS and 
SPIRE instruments onboard the Herschel Space Observatory a 
representative sample of evolved stars, including 78 asymptotic 
giant branch (AGB) stars and post-AGB stars with the aim of 
studying their circumstellar environments and mass-loss history 
by taking advantage of the good angular resolution offered by 
Herschel (with a point-spread function of 5'.'7 full-width at half 
maximum at 70 fim). 

A large variety of shapes (spherical, elliptical, detached, 
and axisymmetric shells) has been encountered by this pro- 
gram. A similar conclusion is reached for smaller angular 
scales (1-10") by the Plateau de Bure/Pico Veleta CO survey 
( ICastro-Can^izoetaLll2 007 . 201 1| . In addition, Mi d-IR and op- 
tical surveys ( Meixner et al... 19991 lUeta et aDl2000h have shown 
that the shells exhibit at least axisymmetry by the end of the 
AGB phase. AGB stars with detached shells in the MESS sample 
have already been discussed by Kerschbaum et al. (2010, 201 j]). 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



and a synoptic paper is in preparation (Cox et al.). The present 
paper focuses on two spectacular cases, namely the O-rich star 
X Her and the C-rich star TX Psc, which show axisymmetric 
shells. 

Extended asymmetric structures around AGB stars have 
been seen in different wavelength bands (which thus probe 
different reg ions of the c ir cumstellar environment): in the 
IR by Youn g et all (119931) lUeta et al .l (2006, 2008, 2010b. 
Ueta (2008), lGeise etal.1 (l2010h. iLadialetal. (2010 ), an d 
Izumiura et al? (!201ll), in H I at 21 cm by Gardan et al. (200^ 
Matthews & Reid (2007), Matthews et al. (2008), Libert etd] 
(2010), and Matthews et al. (2011), in the CO radio l i nes by 
Heske et al. ( 1989) . Kahane & Jura (1996) . iLibert et alJ (l2010h 
and Castro-Carriz o et al.l (1201 ll) . and in the UV bv lMartin et al.l 
(i2007i) . These asymmetries can be triggered by several causes: 

- Wind ejection in a binary system will shape a spi- 
ral stream, according_^ t o hy drodynamical simula- 

tions dl heuns & Joris seij 119931; ^heuns et al. 199J; 

iMastrodemos & MorrisI 119991: ^[Huggins et al) l2009l) . 
The clearest illustration thereof has b een observed in the 
proto-planetary nebula AFGL 3068 dMauron & HugginsI 
l2006l:lMorrisetal.ll2006h . 
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Asymmetric mass loss in single stars, as predicted by dust- 
driven wind models for red giants with cool spots (Reimerd 
2005) and 3D models of carbon stars ( Frevtag & H6fner| 
20081) . may produce an ofF-centred nebula. Such an asym- 



metric mass loss may also impart a kick on the forming 
white dwarf, and evidence th ere of ha s been reported in 
globular clusters (" Hevll l2007allbl 120081: iDavis et al.1 120081: 



Fregeau et al. 2009). 
Matt et alTdlOOOi) demonstrated that a stellar dipole magnetic 



field can focus an initially isotropic wind toward the equa- 
torial plane, thus producing an equatorial overdensity, but 
ISoked (120051) pointed out that in cases where the magnetic 
field plays a global role in the shaping, a binary companion 
is necessary to maintain the field. 
- Interaction between the stellar wind and the interstellar 
medium (ISM), which is the topic of the present paper. 

For X Her and TX Psc, there are strong indications that the 
interaction of their wind with the ISM shapes their axisymmet- 
ric shell. First evidence of this kind of interaction have been seen 
in planetary nebulae, in the form of axisymmetric morp hologies 
oriented along the space motion of the central star (e.g., iJacobv 
198 1: Revnolds 1985; Ki-a utter et al] 119871: iBorkowski etaL 
1990 : Dgani & Soker 19981 and references therein). Among 
non-AGB stars, far-IR parabolic arcs attributable to bow 
sh ocks have been dis covered around a Ori from IRAS data 
by LStencel et al. (1988), and later confirmed by AKARI data 
dUeta et al.. ,20081 Runaway OB stars are a nother c lass of stars 
around which bow shocks are freque nt ( Gull & S ofia 197^ 
van Bur en & McCravl [T988t Ivan BurenI Il993t Ivan Buren et all 
1995). 



fGardan et all ( l2006h . lUeta et al.1 (l2006h . and iMartin et all 
(I2007h mentioned the possibility of wind - ISM interaction 
in relation to AGB mass loss, although the wind velocity is 
very small (< 20 km s"') in comparison with PNe or OB stars 
(~ 1000 km s"') and shock structures seemed therefore unlikely. 

In this scenario, the AGB wind is slowed down as it sweeps 
up ISM material. The swept-up material forms a density en- 
hancement that continues to expand because of the internal pres- 
sure. The thermal emission of the dust at the shock interface 
between t he stellar wind and the ISM has been detected in 
the far IR (lUetaet al.l 120061 l2009l) . or directly in the X/UV, as 
for the bow shock associated with the wind of o Ceti and de- 
tected by GALEX (Martin et al. 2007). Images obtained with 
Spitzer and AKARI have confirmed such an interaction sce- 
nario for th e targets R Hya, R Cas, and o C et (iUeta et al.l 20061: 
lUetalllooi IUetaetalJl2010l) . In RX Lep, iLibertet al.1 (120081) 
found an H I tail extending 0.5 pc southward, as already sug- 
gested for M ira, RS Cnc, and other sources detected in H I 
(Gardan et al. 2006; Matthews & Reid 2007; Martin et aT]|2()07t 
[Matthew s et al.-2008: Libert et al..2010.) . .Sahai & ChronopoulosI 
(2010) presented ultraviolet GALEX images of IRC-i-10 216 
(CW Leo) revealing for the first time its bow shock, which 
is also visible in H erschel PACS and SPIRE infrared images 
(lLadialetalJl2010l) . 



2. Basic data of the targets 

2.1.x Her 

X Her is an oxygen-rich M8III semi-regular variable (SRV) 
of period 95 d an d with a long secondary period of 746 d 
dHinkle et al.ll2002l) . From a radial- velocity monitoring of the 
CO lines at 1.6 /urn, the same authors discovered a period of 
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Fig. 1. The location of X H er and TX Psc in the IRAS col our- 
colour diagram, along with Ivan der Veen & Habind (fT988l) re- 
gions. 



660 d. That all six but one among the SRV stars monitored by 
Hinkle et al. have similar eccentricities and longitudes of pe- 
riastron casts doubts on the orbital origin of the observed ve- 
locity variations. The IRAS Low Resolution Spectrum of X He r 
was classified as E (i.e., silicate emission) bv lKwok et alj(ll997h . 
and as SE6t (i.e., showing a 13 jum feature, most probably due 
to Al -rich oxide dust, and strong amorphous silicate bands) by 
ISioanet al.1 fl996il2003l) . Its position in region Ilia of the IRAS 
colour-colour diagram shown in Fig.[T]is consistent with the LRS 
class referring to a thick dust shell. From a model fit to the IRAS 
data, lYoung et al.l (11993 ) found X Her to be extended to a ra- 
dius of 6.2' at 60 fim; however, the published data (their Fig. 1 1) 
revealed a much narrower profile, with a radius of the order of 
2', and some substructure, which could possibly be attributed to 
the nearby pair of background galaxies UGC 10156a and b (see 
Fig.lHbelow). 

2.2. TX Psc 

TX Psc is one of the brightest and nearest optical carbon stars 
and has been extensively observed from the visual to mm- 
waves. The General Catalogue of Variable Stars (GCVS) indi- 
cates a va riability type Lb , but an average period of 224 d was 
found by Wasatonic ( 1995, 1997), indicating that SRa/b may be 
more appropriate. The photometric variability in the V band is 
of moderate amplitude (0.4 mag peak to peak in the V band; 
see Fig. |2|. However, ten new 2009-2010 observations with 
the HERMES/MERC ATOR spectrograph (Raskin et al. 201^ 
spread over 450 days show large and regular variations in 
the radial velocity from 8 to 15 km s"' (Fig. |2|; these ra- 
dial velocities were derived by correlating the observed spec- 
trum with a template mask containing 3100 lines covering the 
spe ctral range 476 - 654 nm. Earlier velocity measurements 
by iBarnbaum & Hinkld (1 19951) confirm this amplitude range, 
but are too sparsely sampled to reveal any underlying regular- 
ity. Simultaneous AAVSO photoelectric measurements confirm 
the photometric period of about 220 d. Given that the regu- 
lar velocity variations occur on a period substantially longer 
than the SRa/b period of 224 d, these velocity variations can- 
not be interpreted as a manifestation of a pulsating atmosphere 
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Fig. 2. Radial velocities of TX Psc obtained with the HERMES 
spectrograph and an Arcturus cross-correlation mask (black 
squares and left scale), along with AAVSO V magnitudes (open 
triangles and right scale). The horizontal dashed line is the he- 
liocentric centre-of-mass (CoM) velo city from the CO 7=1-0 
radio line (iBarnbaum & HinklelfT99l and references therein). 




Fig. 3. The location of X Her and TX Ps c in the UV d i agram , 
with data for K and M giants taken from iFamaev et al] ( l2005h . 
Black triangles denote stars belonging to the background ve- 
locity ellipsoid, (cyan) open squares to high-velocity stars, and 
(green) open circles to Hercules-stream stars. 



Table 1. Basic kinematical data for X Her and TX Psc . 
Parallaxes and proper motions are from Ivan Leeuwe ^ (I2007h . 
The column labelled 'LSR' provides the kinematical data cor- 
rected from the solar motion. 





XHer 




TX Psc 




Spec. Type 


M8 III 




C5-7,2 (N) 




Var. type'' 


(95 d) 




SR (224 d) 




in 


74.5 




93.3 






47.8 




-55.6 




TU ±cr^ (mas) 


7.30 ± 0.40 


3.63 ± 0.39 






helioc. 


LSR 


helioc. 


LSR 


yu„ cos 6 (mas/yr) 


-68.5 ±0.4 


-57.3 


-33.7 ±0.4 - 


46.4 


fis (mas/yr) 


64.6 ± 0.4 


55.1 


-24.5 ± 0.3 - 


19.0 


RV (km/s) 


-88.95 


-73.6 


13.0" 


13.4 


Syst. vel. (km/s) 


-90'' 




9.3* 




Space vel. (km/s) 


108.0 ±2.0 


89.9 


55.2 ±4.4 


66.9 


Pos. angle (°) 


313.3 


309.0 


234.0 242.6 


Inclination (°) 




-55.4 




13.6 



References: 

a: Average HERMES velocity 

b: CO y = 1 - heliocentric velocity from lBambaum & Hinkld ( Il995h 
and references therein. 



ci pam aev et al. ( 2 00. 
di lKSianefe Jural ( 1199 
e: GCVS 



30l) 



[see iLebzelter et alj ( 1200 Ih for examples of a similar situation 
for radial velocities derived from 4 //m lines]. They could ei- 
ther be an indication of the possible binary nature of TX Psc, 
or yet another example of the so-called long-secondary periods 
encountered in several long-period variables (including X Her; 
[Hinkle et al. 2002; Nicholls et al. 2009) whose origin is still un- 
clear. The radial-velocity monitoring is being pursued to clarify 
the origin of these variations. In the absence of any firm indi- 



cation of the binary nature of TX Psc, we do not consider that 
possibility any further in this paper, especially since it may just 
be of a similar (non-orbital) nature as that observed for X Her 
and discussed above. 

Although the IRAS low resolution spectrum of TX Psc has 
been classified as S (i.e., featureless, stellar-like) bv' Kwok et al.l 
U99 7). its 60 /vm flux is clearly indicative of cool dust, TX Psc 
falUng in region Via of the IRAS colour-colour diagr am (Fig.[l), 
which is the locus of stars with detached dust shells. Young et al.l 
(Il993h indeed found TX Psc to be extended with a radius of 
roughly 3' at 60 /im. From a fit to the near- IR and IRAS spectral 
energy distribution, lEriksson et al.l (Il986h have derive d a dust 
temp erature of 125 K, but a more elaborate study bv IWirsichI 
( 1199 1) suggests that the dust temperature decreases from 1155 K 
at 2.25 steUar radii to 410 K at 30 stellar radii. 



2.3. Kinematical data 

Proper motions, radial velocities, and parallaxes for X Her and 
TX Psc are listed in Table [T] For all calculation s, we used the 
revised Hipparcos data by Ivan Leeuwe B (l200l . and the sys- 
temic velocity derived from radio observations, which is consis- 
tent with the radial velocity derived from optical lines (TablelTJ. 
The space velocity is computed as follows. First the components 
(t/, y, W) of the spatial heliocentric velocity in Galactic coor- 
dinates are computed following Johnson & Soderblom (1987!), 
using the proper motion, the distance, and the radial veloc- 
ity. The sources for these quantities used to compute the space 
velocity are given in Table [Tj The solar motion {U, V, W)q - 
(11.10, 12.24, 7. 25) km s ' ( Schonrich et al.l2010 l) is then added, 
and the a, 6 components of the space velocity, thus corrected 
from the solar motion, are then easily calculated by reversing 
the procedure. In Table [1] the column labelled 'LSR' provides 
all quantities corrected for the solar motion. 

The location of X Her and TX Psc and in the UV diagram is 
shown in Fig. [3] which reveals th at X Her belongs to the group 
flagged as 'high- velocity stars' bv lFamaev et al.l(l2005l) . TX Psc, 
in contrast, happens to fall in the background ellipsoid. 
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Fig. 4. PACS images of X Her. The upper image is at 70 yum, the 
lower at 160 /im. The arrows indicate the direction of the proper 
motion and the position 2000 yrs in the future. The two bright 
dots on the lower left side are a pair of galaxies (UGC 10156a 
and b). The faint source west of X Her is the galaxy MCG+08- 
29-036. 



3. Herschel data analysis 

3. 7. Observations and data reduction 

The images of X Her (Fig. |4]i and TX Psc (Fig. |5]l that we 
present here were obtained by the Photodetector Array Camera 
and Spectrometer (PACS) of Herschel on 2009 December 20 and 
21, respectively, and are part of the Mass-loss of Evolved StarS 
(MESS) guaranteed time ke y programme. Detail s of the MESS 
programme can be found in Groenewege n et aP (1201 ih . which 
describe as well the basic data processing and reduction with the 
Herschel Interactive Processing Environment (HIPE) that was 
used for producing Figs. |4] and |5] All images were oversampled 
by a factor of 3.2, which results in a sampling of 1 pixel per arc- 
sec at 70 fim and 1 pixel per 2 arcsec at 160 jum, to be compared 
with the FWHM of 5'.'7 and 1 1'.'4 at those wavelengths, respec- 
tively. Fig.|6]presents deconvolved images of X Her and TX Psc. 
More details of the deconvolution methods used to produce these 
images can be found in Ottensamer et al. (201 1) . 

3.2. Structure and dynamical age of tine circumstellar nebula 

The images of X Her (Figs.|4]and|6]l and TX Psc (Figs.|5]and|6]l 
reveal a wealth of features that we now discuss. 

For both X Her and TX Psc, there is a noteworthy clump 
in the direction of the space motion (Fig. |6] and bottom right 
panel of Fig. [TOl i. For TX Psc, there is also a ring around the 
star, which is clearly visible in the across scan cut (upper panel 
of Fig. [TOb . This ring is not an artefact of the deconvolution 



Pos. 2000 yrs 



0.16 pc 
2' , 



Fig. 5. Same as Fig.|4]for TX Psc. The bright source east of the 
star is an uncatalogued galaxy. 



process, since it is also visible in the non-deconvolved image 
(Fig. |5]l. Moreover, a similar circular extended (filled) structure 
is also visible in the 70 yum image of X Her (Fig.|6l). TX Psc is 
close to the centre of this ring (the upstream radius is the same as 
the downstream radius; see Fig. [TOl i. The circular symmetry of 
this feature is a clear indication that the mass loss from which it 
originates was isotropic, at variance with the situation currently 
prevailing closer to the star, where a bipolar outflow is observed 
(see Sect.|4|. Moreover, the prominent nature of this ring on the 
PACS images may indicate that it actually corresponds to the 
termination shock, where the reverse shock from the wind-ISM 
interface is facing the freely expanding stellar wind, although the 
available data offer so far no additional support of this hypothe- 
sis. 

The dynamical ages of these features may be estimated 
as follows. The inner ring is at a distance of 17" from 
TX Psc downstream, which corresponds to a kinematic age of 
about 2100 y rs (Fig. [TOb for an a ssumed wind velocity v„, = 
10.5 km s ' (l01ofssoiietalJll993h . neglecting any possible in- 
clination on the plane of the sky. The blob in the direction of the 
space motion is at a distance of 29", corresponding to a kine- 
matic age of 3610 + 410 yrs. The wind-ISM interface limits the 
possibilityof looking further back into the mass-loss history. 

The same procedure and assumptions were applied to the 
heart-shaped X Her With a wind velocity v„ - 6.5 km s ' 
(Gonzalez-Delgado et al. 2003), the bow shock at a distance of 
28" from the star has a minimum kinematic age of 2800±180yrs 
and the end of the tail at 1 12" downstream has a (minimum) age 
of 11200 + 680 yrs (Fig.©. 

Another structural feature in TX Psc deserves discussion. A 
compact feature (representing 2% of the flux in the K band and 
located upstream, at a position angle of 241°) was dete cted by 
the COME-ON -I- observations of lCruzalebes et all (Il998h . but at 
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Fig. 6. Deconvolved PACS images of X Her (left) and TX Psc (right), both in the blue channel at 70 //m. The arrows indicate the 
direction of the proper motion and the position 1000 or 2000 yrs in the future. North is up and east is to the left. 



Fig. 7. The gas density in case E simulation of IWareing et alj 
( l2007ab (corresponding to v»jsm = 125 km s"', rin — 2 cm"^, 
and a mass-loss rate of 5 10"^ yr"'), which reveal the for- 
mation of vortices in the wake of the AGB star. This stage cor- 
responds to the end of the AGB phase, and resembles the X Her 
nebula (compare with Fig.|6l remembering that the viewing an- 
gle is different, since the space motion of X Her is inclined 
by 55° with respect to the plane of the sky). Reprinted from 
MNRAS, Volume 382, Wareing et al., page 382, 2007, with per- 
mission from John Wiley & Sons. 



a distance of only 0'.'35 from the primary source, thus well inside 
the ring described above. Although, at first sight, one would be 
tempted to relate this blob observed close to the star with the one 
observed in the PACS 70 jum image discussed above, they cannot 
be caused by the same physical phenomenon (such as the hydro- 
dynamical instabilities discussed later in the present section), be- 
cause the COME-ON-H secondary source is located withit^ the 




' It is in principle possible that this secondary source is actually lo- 
cated outside the circular ring, but seen inside only in projection. This 
would mean that the secondary source is then located almost along the 
line of sight, at a distance larger than about 4700 AU (= 17"/0'.'0036) 



Fig. 8. The gas density in the Mach 10 simulation of 
Blondin & Koerwer ( 1998). Note the strong blob upstream and 
the highly unstable and inhomogeneous bow shock. This case 
resembles the TX Psc nebula (compare with Fig.|6] remember- 
ing that the viewing angle is slightly different, since the space 
motion of TX Psc is inclined by 14° with respect to the plane 
of the sky). Reprinted from New Astronomy Volume 3, John M. 
Blondin andJoelF. Koerwer, pages 571-582, 1998, with permis- 
sion from Elsevier. 



circular ring observed around TX Psc at a distance of about 17". 
As discussed earlier, structures within this ring cannot be caused 
by the interaction with the ISM. Nevertheless, the inner clump 
is observed at the same position angle as that corresponding to 
the space motion. Another puzzle, possibly related to the for- 
mer, is the correlation existing between the space motion and the 
orientation of the bipolar outflows observed around both X Her 
(where it is almost exactly perpendicular to the space motion) 



from the star. In the absence of any other argument supporting this view, 
this possibility will not be considered further in this paper. 
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and TX Psc (where it seems to be parallel with the space mo- 
tion). The properties of the bipolar outflow are discussed further 
in Sect. |4] Finally, it is rather surprising that the morphology of 
the far-IR nebula does not bear any signature of the inner bipolar 
outflow, especially in the case of X Her where it is perpendicular 
to t he space motion. 

IComeron & Kaped ([l9p^ illustrated the wide variety of 
bow-shock structures that may arise from the interaction of the 
stellar wind with the ISM. They stressed the importance of the 
wind parameters (density and velocity), the star's velocity, and 
the ISM density. These parameters in turn fix the cooling efli- 
ciency of the gas, which appears to be an essential parameter 
in fix ing the struct ure of the wind - ISM interfa ce (.Soker et al.. 
[19911 Iva n Buren 1 991 IComeron & KaDei]ll998h : if the cooling 
is eflicienj^ the gas can be compressed to high densities, and 
the shocked layer remains very thiij^- This situation is described 
well by the isothermal approximation. On the contrary, if cool- 
ing is inefficient, the high temperature of the gas limits the com- 
pression factor, which for the li miting case of no cooling , would 
be equal to 4 (abiabatic shock; Comeron & Kaper'1998'). In this 
situation, the shocked layer is thick. Different instabilities are en- 
countered in thin and thick shocks (see lComeron & Kapeilfr998l 
for a review). 

In thin shells, iDeani et alJ (Il996bllah studied the so-called 
transverse acceleration instability, which occurs as a conse- 
quence of the acceleration in the flow normal to the shock caused 
by its curved surface. Another instability appearing under the 
same conditions is the non-linear thin-shell instability, which 
arises from the shear in shock-bounded slabs produced by devia - 
tions from equilibrium ( Vishniac 1994; Blondin & Marks 1996). 
If there is a shear along the surface of the bow shock, Kelvin- 
Helmholtz instabilities will develop as well. Numerical simu- 
lations relative to bow shocks with thin shells were performed 

by_ 

, and lBlondin & Koerweri (119981) . 
If the shocked material cools rather inefficiently, a thick shell 
will form, which will be subject to other types of instabilities, 
mostly the Rayleigh-Taylor (R-T) instability, because the den- 
sities of the gas entering from either side of the shock will not 
match in the absence of cooling. When two fluids of different 
densities are penetrating each other, R-T instabilities manifest 
themselves in the form of mushroom- and finger- shaped pro- 
tubera nces of the lighter fluid into the heavier one. ISoker et all 
demonstrated that R-T instabilities occur when a stellar 
wind or planetary nebula shell is colliding with a low-density 
ISM with a high relative velocity (v,jsm > 100 km s"'), leading 
to a fragmentation of the shell. A bump in the upstream direc- 
tion, attributed to a R-T instability developing at an early stage 
of the wind - ISM in teraction, is seen fo r instance in the plane- 
tary nebula NGC 40 (iMartin et al.ll2002l) . At a later stage, as the 
shell fragments, Kelvin-Helmholtz (K-H) instabilities appear as 



byHo ker et al. (1991), Dgani et al. ( 1996b), Comeron & Kaper 
([l998h .; 



^ See Soker et al. (1991) for a basic description of the relation be- 
tween the wind/ISM properties and the cooling efficiency. Dense gas 
generally gives rise to efficient cooling. 

' Wilkin fitting - geometrical fittin g of th e bow shock shape based on 
the analytical formulation by WilkinI dl996h . which will be attempted in 
Sect. l3.4l - applies to the thin-shell approximation. 

'* A related i nstability is the Richtmyer-Meshkov instability (e.g., 
lNeedharnll201(]t) . which develops when a shock wave hits the interface 
of two fluids with different densities, such as when the new shell of 
an expanding wind bubble hits the already established bow shock. This 
may happen as AGB stars lose mass in a succession of discrete events, 
as suggested by th e successive shells o bserved around e.g. the carbon 
Mira IRC +10216 dMauron & Hugginsll200Ql . 



a sequence of vortices peeling off the bow shock and moving 
downstream (I Wareing et al.ll20d7allbl) . 

The nebulae around X Her and TX Psc show very un- 
usual features, such as the upstream blob and clumps visible 
along what seems to be a bow shock. Investigations of the 
shape of the bo w shock have been made us i ng hydrodynamical 
simulations by Blondin & Koerwer (1998), Comeron & Kape^ 
(1998), Viflaveretal. (2003), and Wareing et al (2007a). The 
Case E simulation of the last set of authors (corresponding to 
v*,isM - 125 km s"', nn - 2 cvoT^, and a mass-loss rate of 

5 X 10"^ M0 yr"^; see Fig.|7]i, at stage a (corresponding to the 
end of the AGB phase), looks very similaf] to the X Her nebula 
(compare Figs. |6] and |7]), since it combines an upstream shock 
surface bent towards the star, and K-H vortices in the down- 
stream wake. Interestingly, the 90 km s ' LSR space velocity 
of X Her (Table [T]) is similar to the velocity used in the Wareing 
et al. simulation, but the current mass-loss rate is a factor of at 
least 30 lower than the Wareing value (see Sect. 14. 11 1. On the 
other hand, the TX Psc nebula bears similarities (especially the 
strong upstream blob, and the inhomogeneous bow) with the 
large Mach number (10) and isothermal simulations of Blondin 

6 Koerwer (1998; compare Figs. |6] and [8]). A large Mach num- 
ber is indeed expected in the present situation, since the X Her 
space velocity is 67 km s ' and the sound speed is of the order of 
1 km s ' (for a temperature of the order of 80 K, as inferred from 
the analysis of Sect. 13.31 ). The nature of the ring observed in the 
TX Psc image is unclear Similar rings are clearly visible in the 
Wareing et al. simulations when the mass-losing star has reached 
the post-AGB phase (thus posterior to the stage displayed in 
Fig.|7]). The ring thus corresponds to the planetary nebula shell. 
This cannot be the case for TX Psc, although the ring possibly 
traces a discrete episode of strong mass loss. Alternatively, it 
could be the termination shock. 



3.3. Interpreting the emission of ttie upstream blob 

As noticed at the end of Sect. 13.21 an upstream blob is 
clearly visible in the Mach- 10 isothermal-shock simulations of 
Blondin & Koerwer (.1998.) (Fig.lSt. as w ell as in several of the 
simulations of Comeron & Kapeir(ll998h . Such a blob is present 
in the nebulae associated with both X Her and TX Psc. To inves- 
tigate the nature of the emission associated with this blob, we 
performed aperture photometry in the 70 and 160 fim bands, for 
TX Psc, where it is cleary detached from the surrounding ma- 
terial (lower right panel of Fig. [TOb . A radius of 5 arcsec was 
adopted to fully cover the blob. The background flux was esti- 
mated, from the periphery of the image, to be 3x10""* Jy per pixel 
at 70 jum and 5 x 10""^ Jy per pixel at 160 fim. With these values, 
we get total fluxes of 0. 191 Jy at lO/jm and 0.039 Jy at 160 yum, 
or a ratio Fv,7o/^v,i60 = 4.89 (see also the upper panel of Fig.fTTl 



' We stress, however, that Wareing et al. numerical models do not in- 
corporate dust grains, and only handle gas warmer than 10"* K (because 
of the lack of an adequate cooling function for gas with lower temper- 
atures). In such environments, dust grains would not be able to survive 
and therefore far-IR dust thermal emission is not expected. Although 
Fig. |7] shows structures resembling what is seen in PACS images, this 
model traces the hot gas density, which is not expected to emit the ob- 
served far-IR emission. Despite these shortcomings, the striking sim- 
ilarity beween the PACS and model images could be accounted for if 
the far-IR radiation comes from atomic emission lines in shocked gas. 
Alternatively, one may speculate that future models including the proper 
gas cooling function will find similar structures involving gas at lower 
temperatures, with dust naturally forming in regions of high gas densi- 
ties and moderate temperatures. 
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Fig. 9. Intensity cuts for X Her from the PACS 70 //m decon- 
volved image (right panel of Fig. O along two perpendicular di- 
rections through the star, one corresponding to the space motion 
direction (lower panels). The dynamical ages on the upper scale 
correspond to a wind expansion velocity of 6.5 km s 
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Fig. 10. Same as Fig. [10] for TX Psc and a wind expansion ve- 
locity of 10.5 km s"'. 



which compares the intensity profiles in the 70 and 160 fim 
bands). The observed light in this blob could in principle come 
from emission lines (such as [O 1163 fim, [O 1] 145 fim, and [C 
II] 158 /im) excited by the shocly. However, Spitzer/IRS spec- 
troscopic observations of the R Hya bow shock did not detect 
atomic lines (Ueta,201 1) . and the detected far-IR radiation from 
the bow suggests instead the presence of very cold (20-40 K) 
material (Ueta et al."2009^. The same conclusion was reached 
by van Buren & McCray ( 1988) and Dgani et al. ( 1996a) in the 
context of OB runaway stars; they note that the dust trapped in 
the bow shock shell reradiates efficiently about 1% of the star's 
bolometric luminosity at far- infrared wavelengths. 

Therefore, we make here the working hypothesis that the 
light emitted by the blob is from dust emissiorQ, and derive the 
corresponding dust temperature. From Kirchhoff 's law, we may 
write the dust emissivity jy as 



jv = ky JbbAT), 



(1) 



^ Remember that H„ emission as well as UV radiation (caused by 
excited H2 molecules) are someti mes associated with AGB wind - ISM 
interaction, as i n the case of Mira jMartin et alj|2007l) and R Hya jUetal 
l2008l :IUetaetal. 2008) 

^ Solving this issue would require the acquisition of a spectrum for 
this blob. 
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Fig. 11. Intensity plots for X Her and TX Psc. Upper plot: in- 
tegrated intensity of X Her in the direction of the northern blob 
(PA: -7° to -17°) from the image of Fig. H) Lower plot: in- 
tegrated intensity for TX Psc in direction of the blob/proper 
motion (PA: 228° to 250°), from the image of Fig. |5] In both 
plots the (blue) continuous curve corresponds to 70 fim, the (red) 
dashed curve to 160 //m. The distances used to convert the angu- 
lar scales into distances in AU are listed in Table [T] 

where JbbAT) is the black body emissivity at frequency v and ky 
is the dust absorbing coefficient. Since TX Psc is a carbon star, 
we adop ted the absorp tion coefficient for amorphous carbon as 
given bV' Drain3(ll98lh . who provides the Q factors (ratio of the 
absorption cross section to the geometrical cross section) as a 
function of wavelength. Since ky is proportional to Q and we are 
only interested in ratios, we may write 



v,70 



',160 



Gto Jbb,vO0)(T) 

2l6() yBB,v(160)(r) 



where 270/2160 = 512/219 
is that satisfying the relation 



(2) 



2.34. Thus the dust temperature 



: 4.89/2.34 = 2.09. 



(3) 



7BB,v(70)(r) 
jBB,v(i60)iT) 

A temperature of 8 1 K for the blob satisfies the above condition. 
Since the absorption coefficient is for dust containing only amor- 
phous carbon, the temperature derived above must be seen as a 
rough estimate. The same analysis for X Her, adopting the ra- 
tio Q70/216O - 185/25 - 7.4 for astronomical silicates (iDraind 
1985), yields a temperature of 40 K. The composite colour im- 
age of X Her presented in Fig.[T2]does not show evidence of any 
large temperature variations in the nebula, which is dominated 
by the emission at 70 jum. 



3.4. Wilkin fitting of tine bow shock interface 

The far-IR parabolic-like arc structures around X Her and 
TX Psc can be interpreted as the interface regions between 
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Fig. 12. Composite image merging the PACS 70 and 160 fim 
images of X Her (respectively in the green and red channels). 



the ISM and AG B winds , as already rec ognized for R Hya 
(|Uetaetal.ll2006h. R Cas dUeta et al.ll20"Toh . and IRC +10 216 
(iLad ial et al.'2010). Assuming that (1) the ram-pressure balance 



is established between the ambient ISM and AGB winds and (2) 
momentum is conserved across a physically thin shock interface 
(i.e. the shock is radiatively cooled efficiently), one can express 
the bow shock shape analytically as a function of the latitudinal 
angle measured from the apex of the bow with respect to the 
position of the central star as formulated bv .Wilkin.(,1996i) 



R(e) = R^ 



V3(l -0 cot 61) 



sinO 



(4) 



where Rq is the stand-off distance between the star and bow apex, 
and defined as 



^0 = 



\4;rpj 



(5) 



ISM 



for which M is the mass-loss rate, is the isotropic stellar wind 
velocity, pisM is the ambient ISM mass density, and v, jsm is the 
relative velocity of the star with respect to the ambient ISM. 

Because the observed surface brightness enhancement rep- 
resents approximately a conic section of an inclined bow shock 
paraboloid by the plane of the sky intersecting the central starj^ 
the inclination angle and deprojected stand-off distance of the 
bow shock paraboloid can (in principle) be determined by fit- 
ting the apparent shape of the bow shock using the analytical 
Wilkin function (the detailed description of the Wilkin fitting 
procedures may be found in Cox et al., in preparation). 

As discussed in Sect. 13.21 the Herschel images of the targets 
at higher spatial resolution have resolved structures in the bow 
shock interface that resemble vortices caused by in stabilities at 
the contact discontinuity (e.g..^ Wareing et afl2007bh . In particu- 
lar, the X Her image shows that the bow shock interface is bent 
in the direction of the star (Fig.|6]). This bent-interface structure 
appears typical of bow shocks caused by a large relative veloc- 
ity between the star and the a mbient ISM, v.i ,. lsm, a s shown in 
hydrodynamical simulations bv lWareing et al.l (l2007 b) (Fig.|7]i. 



^ Assuming the observed surface brightness is mainly due to ther- 
mal dust emission, the surface brightness is roughly proportional to the 
dust column density along the line of sight in the optically thin limit, 
which is true for the far-IR wavelengths. If the dust distribution is an 
inclined paraboloid, the column density becomes the greatest at the in- 
tersection of the structure with the plane of the sky and therefore the 
observed surface brightness distribution represents the conic section of 
the paraboloid. 



The presence of these bent structures obviously makes the 
application of the Wilkin fitting problematic. Nevertheless, we 
performed the Wilkin fitting of the bow shock with and without 
the bent region and obtained essentially the same results. While 
we present the results of Wilkin fitting below assuming that a 
conic section of the bow shock paraboloid is still recoverable 
from the overall shape of the observed far-IR surface brightness 
distributions, we recall that the observed bow shocks of X Her 
and TX Psc show asymmetric structures that were not seen in 
any previous AGB wind-ISM interaction cases. 

Using the method outlined by Cox et al. (in preparation), the 
Wilkin fitting then yields a best-fit bow shock paraboloid with 
Ro = 36" (corresponding to 4930 AU at the parallactic distance 
of 137 pc) at the position angle (PA; defined as degrees E of N) of 
330° for X Her and with Rq = 37" (corresponding to 10175 AU 
at the parallactic distance of 275 pc) and at the PA of 238° for 
TX Psc. 

The PA obtained above (characterizing the relative motion 
of the AGB star with respect to the ambient ISM as determined 
from the apparent heliocentric orientation of the bow shock 
paraboloid, and denoted v*.ism) must be compared with the PA of 
the stellar space motion in the LSR (v» lsr), as listed in Table [T] 
They are, for both stars, roughly consistent with each other: 330° 
and 309° for X Her, and 238° and 243° for TX Psc. The 20° dis- 
crepancy observed for X Her may be an indication that the ISM 
around that star is not at rest in the LSR; however, we do not 
attempt to infer the velocity of that ISM flow, because of the un- 
certainties associated with Wilkin fitting for such peculiar shock 
shapes as observed around the two stars under consideration. 

Assuming as a first approximation that the ISM is at rest 
in the LSR, (thus v» ism - v,,lsr, the latter value being listed 
in Table [TJ, Eq. |5] may then be used to derive the ISM density 
close to X Her and TX Psc, knowing their mass loss rate and 
wind velocity. Adopting for the se values M = 1.5 X 10"^ Mp 
yr ' and 6.5 km s ' for X Her (iGonzalez-Delgado et aLll2003l) 
and M ^ 9.1 x 10' *^ Mq yr^' and 10.5 km s"' for TX Psc 
(lOlofsson etal.lll993h . we find = 0.67 and 0.16 atoms cm"^ 
for X Her and TX Psc, respectively. These values appear quite 
reasonable, given the high Galactic latitudes of these two stars 
(Table©. 



4. Bipolar outflows and global picture 

4.1. XHer 

An extensive descripti on of the cir cumstellar environ ment of 
X Her can be found in lGardan et al.l (^006) and Matthew s et al.l 
( 1201 Ih . and only the most noteworthy features will be 
repeated here. X Her has distinctly double-component 
lines in CO {J -2 - 1) and (7 = 3 - 2), i.e., a narrow 
feature centred on a broader feature, w ith Vw ~ 3 and 
~ 9 km s ', respectively ( Kahane & Jur^ [T996^ Knapp et al 



1998 
2002; 



Kerschbaum & OlofssonI | 1 , 999| ; Olofsson et al 



20031: iNakashimai ,20051 



G onzalez- Delgado et al.i . _ 

iCastro-Ca rrizo et al.l l201 ih . iGonzalez-Delgado et all (l2003h ob- 
serve a similarly double-component SiO(v = 0, 7 = 2 - 1) line; 
from CO data, they derive mass-loss rates of 1.5 x 10"^ Mq y"' 
and 0.4 x 10"^ M© y"', respectively, for the broad and narrow 
components. Although iKnapp et al.l (1 19981) argued that this 
complex line profile may be caused by episodic mass loss with 
highly varying gas expansion velocities, it seems more likely 
now that it results from a complicated flow geometry. 

The s patial information prov ided by the CO radio line obser- 
vations of iKahane & Jural d 19961) suggests that the broad plateau 



Jorissen et al.: Herschel observations of X Her and TX Psc 



9 



is a bipolar outflow oriented in the SW - NE direction, with 
the SW side facing towards us, a structure that was confirmed 
by Nakashima (2005). The most recent CO radio fine interfer- 
ometry data of ICastro-Carrizo et alJ (|201 ih clearly resolves an 
hour-glass structure in the higher-velocity emission with a po- 
sition angle of 45°, i.e., this indicates that a collimated bipolar 
wind in this direction has shaped the original AGB circumstel- 
lar envel ope. The inclination angle of the outflow is likely to 
be small. ICastro-Carrizo et akl (1201 ll) find no evidence of a (ro- 
tating) disk perpendicular to the outflow axis. We obtained SiO 
(v - 0,J = 1-0) line data towards X Her using the Very Large 
Array (details about these data and their reduction are given in 
the Appendix). These data probe a region smaller than that of the 
CO data because the SiO molecules are confined to a smaller re- 
gion and they are eff'ectively excited only in the inner more dense 
parts of the circumstellar envelope. Nevertheless, the channel 
map data presented in Fig. [13] show a structure very similar to 
that of the CO data, a bipolar outflow at position angle 45° and 
a hint of an hour-glass structure. The position-velocity diagram 
along the outflow axis is also very similar to that obtained from 
the CO line data (Fig.fTJb. Interestingly, the position- velocity di- 
agram perpendicular to this axis shows a structure that could be 
interpreted in terms of a rotating disk. 

iGardan et all (l2006l) and lMatthews et all (l201lh detected the 
H I emission at 21 cm from the circumstellar shell around X Her, 
and found it to be extended al ong the direction of the space mo- 
tion (> 6' or about 0.24 pc). iMatthews et alJ (l201lh found, in 
their H I data and on an arcmin scale, some evidence for the con- 
tinuation of the bipolar outflow observed on small scales (+10") 
in the molecular CO and SiO lines. 

The surprising picture that emerges from the combination 
of all the data, including the new Herschel data set, is that of 
a bipolar outflow in the innermost region (+10"), that appears 
to be almost exactly perpendicular to the far-IR and H I nebula 
oriented along the direction of space motion and caused by the 
interaction between the stellar wind and the ISM. It is not at 
all clear whether this is just a chance coincidence, or whether 
there could be a causal relationship between the orientation of 
the bipolar outflow and the space motion. The absence of any 
sign whatsoever of the bipolar outflow visible in the Herschel 
images would tend to dismiss the latter possibility. The existence 
of a bipolar outflow seen in CO and SiO molecular lines raises 
another question: how can this outflow be reconciled with the 
bow shock seen all around the star and signaling that mass loss 
is (was?) roughly isotropic? Two solutions may be envisioned: 
(i) even in the presence of a bipolar outflow, mass continues to 
be ejected in all directions, albeit probably at a lower rate; (ii) 
the bipolar outflow has formed more recently than the features 
seen in the far IR. 



4.2. TX Psc 

The close environment of TX Psc, as probed by the molecular 
lines, appears to be quite complex. The observed line profiles of 
the CO line emission towards TX Psc are indeed peculiar, as first 
noted by Heske et al. ( 1989). As a consequence, the expansion 
velocity derived from the observed ra dio line profiles range from 
7.5 to 12.2 km s'' ( see Table 2 of ISchoier & OlofssonI lIOOll 
and iLoup et aL 1993h de p endin g on the transition and the tele- 
scope used. ' Heske et al.l (1X989) mapped the circumstellar CO 
line emission around this object and interpreted the observed 
asymmetries as being produced by either a bipolar outflow or 
a highly clumped wind. No detailed modelling exists for the cir- 
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Fig. 13. Maps from the Very Large Array in the SiO (v = 0, 7 = 
1-0) line for X Her The upper panels show 2 kms ' wide 
channel maps, and the lower panels show position (in arcsecs, 
along the vertical axis) - velocity (in km/s, along the horizontal 
axis) diagrams along the NE - SW [labelled (+4,+4) to (-4, -4)] 
and NW - SE [labelled (+4,-4) to (-4, +4)] directions. Velocities 
are with respect to the local standard of rest (systemic velocity 
of X Her is -72.5 km/s). The lower panels show evidence of a 
bipolar outflow in the NE - SW direction (thus perpendicular to 
the space motion) and of a rotating disk in the NW - SE direc- 
tion. The contour increments are 10 mJy km s ' per beam for 
the channel maps in the upper panel and 5 mJy per beam for the 
position-velocity diagrams in the lower panel. The beam size is 
0.45" X 0.40" (with a position angle of -88 deg). 



cumstellar shell around TX Psc, and its evolutionary status re- 
mains uncertain. 

The asymmetry detected by iHeske et"^ (Il989t) and possi- 
bly attributable to a bipolar outflow is oriented along the NE - 
SW direction, with the red- and blue-shifted part of the emission 
separated by about 10 arcsec. However, there is some ambiguity 



10 



Jorissen et al.: Herschel observations of X Her and TX Psc 



about this direction in the original paper. Inter estingly enough, 
both lunar occultation data ( Richichi et al. and Come-On+ 
adaptive-optics observations (Cruzalebes et al. 1998) claim as 
well that the source is asymmetric. The September 1994 Come- 
On+ observations in the K band clearly reveal a secondary 
source 0.35" to the SW (with a position angle of 241°) of the 
primary source, with a flux ratio of 2%. The lunar occultation 
observations indicate that the brightness profile of TX Psc has a 
structure that becomes increasingly more complex when going 
from the V band to the L band, where it is double peaked (on an 
angular scale of a few milliarcsecs). The V band reveals some 
secon dary peaks on e ither side of th e central core. In the visual 
band. lHartkopf et all ([l996.) and Mason etal.l (1 1999.) could not, 
however, resolve TX Psc between 1993 and 1997 with a resolv- 
ing limit of 0.054". 

The gas mass-loss rate of T X Psc has been derived by several 
authors. ADlofsson et al.l (Il993h found a very moderate gas mass- 
loss rate of 9. 1 X 10~^ M© yr"^ fro m CO observation s and a wind 
expansion velocity of 10.5 km s'. ILoup et al. (Il993h found5.6x 



lO"'' Mq yr-i and in a later paper ISchoier & OlofssonI (1200 lb 



derived an expansion velocity of 7.5 km s~\ with a complex 
CO line shape that prevented them from attempting to derive the 
mass-loss rate. 

The new Herschel-PACS data described in Sect. [3] confirm 
the presence of an axisymmetry in the far IR with a NE - SW 
axis and on an arcmin scale. However, in contrast to the situation 
prevailing for X Her, the asymmetries observed around TX Psc 
on an arcsec scale have the same orientation as the proper mo- 
tion and the far-lR nebula. Nevertheless, an almost circular shell, 
with a radius of about 17", surrounds TX Psc and does not bear 
any signature of the asymmetries observed on a much smaller 
angular scale. 

5. Conclusions 

The PACS images of the O-rich star X Her and the C-rich star 
TX Psc obtained by Herschel at 70 and 160 fim reveal a complex 
structure with signatures of interaction between the AGB wind 
and the ISM. Their space velocities, corrected from the solar 
motion, are 90 and 67 km s respectively. Its v elocity locates 
X Her among the high-velocity group defined by 'Fama ev et alJ 
(l2005h in their kinematic study of giant stars from the Hipparcos 
Catalogue. The most noteworthy feature is, for both stars, the 
presence of a clump located right along the direction of the space 
motion. It is very likely the signature of instabilities arising in 
the shock front between the ISM and the AGB wind. The kine- 
matical age of these clumps is (at least) about 2000 yrs for both 
stars. Although emission lines excited by the shock may con- 
tribute to the observed infrared light, assuming instead that only 
dust emission is involved yields a dust temperature of about 40 K 
for X Her and 80 K for TX Psc. The definite identification of the 
physical process at the origin of the light emitted by the bright 
clumps must await their detailed spectroscopic analysis. Further 
downstream along the bow shock, the PACS images reveal the 
presence of fainter clumps that may be associated with Kelvin- 
Helmholtz vortices peeling off the bow shock and moving down- 
stream. There is actually a striking similarity between the PACS 
image of X Her, and the predictions of the hydrodynamical sim- 
ulation of1^^ejng|eral]j2^07a) for the case v, ism = 125 kms' 
(as compared to the observed value of 90 km s '), uh -2 cmT^, 
and a mass-loss rate of 5 x 10"^ M© yr"' (a factor of 30 larger 
than the current mass loss rate, however), where a front shock 
is bent towards the mass-losing star along the upstream direc- 
tion, and K-H vortices are downstream along the bow shock. The 



strong blob visible upstream of TX Psc may be an indication that 
the wind and the ISM flow are highly supersonic, according to 
the simulations of Blondin & Koerwer ( 1998). 

The Wilkin fitting of the bow shock shape is made difficult by 
its disturbed, highly structured shape (K-H instabilities?) and its 
reverse bending centred on the upstream direction. Nonetheless, 
the position angle of the bow shock obtained by a Wilkin fitting 
is roughly consistent with the position angle of the space motion 
of the star with respect to the LSR (the largest discrepancy is 20° 
for X Her), indicating that the ISM is not far from being at rest 
in the LSR local to the stars. 

A very puzzling conclusion of the present analysis is that for 
both stars the axisymmetry caused by the space motion and its 
interaction with the ISM seems to be maintained down to very 
small angular scales, even beyond the circular feature present on 
intermediate scales in both TX Psc and X Her: 

- In X Her, a secondary source almost perfectly aligned with 
the space motion was discovered by earlier COME-ON+ ob- 
servations, despite being located well within the ring ob- 
served with PACS. 

- Around TX Psc, a bipolar outflow extending over a few arc- 
seconds has been reported by radio observations (with some 
ambiguity however, about its orientation, but seemingly par- 
allel to the space motion). A bipolar outflow was also re- 
ported for X Her, and for that star, the space motion and the 
outflow are clearly perpendicular to each other. 

The existence of those correlations between the axisymme- 
try present in the inner and outer circumstellar regions of TX Psc 
and X Her (despite circular structures being present on interme- 
diate scales, possibly associated with the termination shock) is 
the most puzzling conclusion from the present work, for which 
we do not currently have any convincing explanation. Similarly, 
it is rather surprising that the far-IR nebula contains a nearly per- 
fectly circular shell (either interior to or identical to the termi- 
nation shock), indicative of (past) isotropic mass loss, whereas 
there is at the same time an inner bipolar outflow. Could this be 
an indication that the mass loss became bipolar quite recently? 
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Appendix A: SiO line VLA observations of X l-ler 

X Her was observed in the rotational lines of SiO using the Very 
Large Array in C configuration. As part of a 10-hour observing 
run on 2001 August 24, about one hour was spent to observe the 
{ J - 1 - 0, V = 0) transition at 43.424 GHz in dual circular po- 
larization. A long run of nine hours in right-hand circular polar- 
ization (RCP) was observed on 2004 March 2004. We combined 
these data with an additional hour dual polarization observations 
taken on 2004 April 22. 

All these observations used a total bandwidth of 6.25 MHz 
(~ 40 km s"') centred on a Doppler velocity of -71 km s"' with 
respect to the Local Standard of Rest. The bandwidth was di- 
vided in 64 spectral channels, corresponding to a channel sepa- 
ration of ~ 0.67 km s ' (97.7 kHz), for all three observation runs. 
Also common to the observations was the fast switching mode 
to the phase reference source J1549+5038 using a typical cycle 
of ~40 and 100 seconds for calibration and target source. Using 
3C286 as a flux calibrator, we measured fluxes for J1549+5038 
of 0.52 (±0.01) and 0.61 (+0.01) Jy for 2001 August and 2004 
March, respectively. We set it to be 0.62 Jy for the 2004 April 
observations. 

The data were reduced using standard editing, calibration, 
and imaging procedures available in the Astronomical Image 
Processing System. The RCP data was combined with the dual 
polarization data to Stokes I, assuming that the thermal radiation 
is unpolarized (i.e.. Stokes V = 0). The concatenated visibility 
data set was imaged into a cleaned cube of ~13" x 13" over 
-87 to -54 km s"' with an angular resolution of about 0.5" and 
a RMS noise of typically 2.2 mJy per channel. The results of 
these observations are displayed in Fig. [13] 



